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Abstract 

The critical behavior of an MiV-component order parameter Ginzburg-Landau 
model with isotropic and cubic interactions describing antiferromagnetic and 
structural phase transitions in certain crystals with complicated ordering is stud- 
ied in the framework of the four-loop renormalization group (RG) approach in 
(4 — 2e)-dimensions. Using dimensional regularization and the minimal sub- 
traction scheme, the perturbative expansions for RG functions are deduced for 
generic M and N and resummed by the Borel transformation combined with a 
conformal mapping. Investigation of the global structure of RG flows for the 
physically significant cases M = 2 and iV = 2, A" = 3 shows that the model has 
a three-dimensionally stable fixed point different from the Bose one. The critical 
dimensionality is proved to be exactly two times smaller than its counterpart in 
the real cubic model: = \N^. The numerical value = 1.447 ± 0.020 is 
obtained from resumming the known five- loop e-series for N^. Since < 2, 
the critical thermodynamics of the model relevant to the phase transitions in real 
substances should be governed by the complex cubic fixed point with a new set 
of critical exponents: 7 = 1.404(25), v = 0.715(10), rj = 0.0343(20) for N = 2 
and 7 = 1.390(25), v = 0.702(10), n = 0.0345(15) for N = 3. 
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We study the critical behavior of an MA^-component field model with a cubic 
anisotropy having a number of interesting applications to phase transitions in three- 
dimensional simple and complicated systems. The Landau- Ginzburg- Wilson Hamilto- 
nian of the model reads: 
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where Lp a , for each a, is the M-component field in D = 4 — 2s dimensions, mo, uq and 
v are the "bare" mass and coupling constants, respectively. 

For M = N = 2 Hamiltonian (JJ) describes the structural phase transition in Nb02 
crystal and the antiferromagnetic phase transitions in TbAu2 and DyC2- Another 
physically important case M = 2, N = 3 is relevant to the antiferromagnetic phase 
transitions in ^ItCIq, and Nd crystals The magnetic and structural phase 

transitions in a cubic crystal are governed by model (JI]) at M — 1 and Af = 3 0. In 
the replica limit A" — > (M = 1) Hamiltonian ([!]) is known to determine the criti- 
cal properties of weakly disordered quenched systems undergoing second-order phase 
transitions with a specific set of critical exponents 0. Finally, the case M = 1 and 
A" — > 00 corresponds to the Ising model with equilibrium magnetic impurities ||. In 
this limit the Ising critical exponents take the Fisher renormalization ||. Since the 
static critical phenomena in a cubic crystal as well as in randomly diluted Ising spin 
systems are well understood |7|, §, || [10], [11], [12], we will focus here on the critical 



behavior of the above mentioned multisublattice antiferromagnets. This is the case of 
M = 2 and N = 2,N = 3 in fluctuation Hamiltonian (g). 

For the first time the magnetic and structural phase transitions in crystals with 
complicated ordering described by model (0) were studied by Mukamel and Krinsky 
within the lowest orders in e . A new three-dimensionally stable fixed point ("unique" 
point), different from the Heisenberg or the Bose one, was predicted. The point was 
shown to determine a new universality class with a specific set of critical exponents. 
However, for the physically important case N = 2, the critical exponents of the unique 
fixed point turned out to be the same as those of the Heisenberg one within the two- 
loop approximation. For the years an alternative analysis of critical behavior of the 
model, RG approach in fixed dimensions, was carried out within the two- and three-loop 
approximations JL4|, T5| . Those investigations gave the same qualitative predictions: the 
unique stable fixed point does exist on the 3D RG flow diagram. However, the critical 
exponents computed at this point with the use of different resummation procedures 
were proved to be close to those of the Bose fixed point rather than the Heisenberg 
one. It was also shown that both unique and Bose fixed points are very close to each 
other on 3D diagram of RG flows, so that they may interchange their stability in the 
next orders of RG approximations ||15|| . Recently critical properties of the model have 
been analyzed in third order in e [|T^, [T7]]. Investigation of the fixed points stability 
and calculation of the critical dimensionality A^, of the order parameter separating 
two different regimes of critical behavior confirmed that the model (jl]) possesses the 
anisotropic stable fixed point for N = 2 and A^ = 3. By making use of a new approach to 
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summation of divergent field-theoretical series JT8J based on the Borel transformation 
in combination with a conformal mapping the critical exponents estimates for 
the unique stable fixed point were obtained [17]. The values appeared to be close to 
those of the Heisenberg point in contradiction to the numerical results given by RG 
approach directly in 3D [TJ], |15| . The cause of such a distinction in estimates of the 
critical exponents is the quite short, three- loop e series used. So, the aim of this work 
is to extend existing three-loop e expansions of the model to the four-loop order and to 
study critical phenomena in substances of interest more carefully. Namely, on the basis 
of the four-loop expansions for RG functions obtained in the framework of dimensional 
regularization and the minimal subtraction scheme we analyze the stability of the fixed 
points and calculate the critical dimensionality of the order parameter field ■ Then 
we give more accurate critical exponents estimates applying the summation approach 
developed in Ref. fl8|] to the four-loop series. 

The four-loop e expansions for the /3-functions of the model are as follows 
4 1 
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where £(3),£(4), and ((5) are the Riemann ( functions. 

From these equations we find the formal series for the fixed points. Then we cal- 
culate the stability matrix eigenvalues using an approach based on the Borel transfor- 
mation 

F(2e;a,b) = ±A k (\) f e"& (— )V— ) (4) 
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modified with a conformal mapping u = ||19|| , which does not require the 



knowledge of the exact asymptotic high-order behavior of the series ||18|| . If both 
eigenvalues are negative, the associated fixed point is infrared stable and the critical 
behavior of experimental systems undergoing second-order transitions is determined 
only by such a stable point. For the most intriguing the Bose and the complex cubic 
fixed points our numerical results are presented in Table I. It is seen that the complex 
cubic fixed point is absolutely stable in D = 3 (2e = 1), while the Bose point appears 
to be of a " saddle- knot" type. However, A2's of either points are very small on the 
four-loop level, thus implying that these points may swap their stability in the next 
order of RG approximation. 

In addition to the eigenvalues, we also calculate the critical dimensionality of 
the order parameter. The four-loop expansion is 

iVf = 2 - (2e) + — [6C(3) - ll (2e) 2 + — |45C(3) + 135C(4) - 600C(5) - ll {2sf. (5) 
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Instead of processing this expression numerically, we establish the exact relation = 
which is independent on the order of approximation used. Here and are 
the critical dimensionalities in the complex and in the real cubic model, respectively. 
The five-loop e-expansion for was obtained in Ref. ||. Resummation of that series 
gave the estimate N* = 2.894(40) 0. Therefore we conclude that iVf = 1.447(20) 
from the five-loops. Practically the same estimate = 1.435(25) follows from a 
constrained analysis of that takes into account that in two dimensions N^- = 2 . 
So, the phase transitions in the NbC>2 crystal and in the antiferromagnets TbAu2, DyC% 
K2lrCle, TbD2, and Nd are of second order and their critical thermodynamics should 
be controlled by the complex cubic fixed point with a specific set of critical exponents, in 
the frame of given approximation. Corresponding four-loop critical exponents estimates 
are displayed in Table II. 

Although our results seem to be self-consistent, still there is definite contradiction 
with the nonperturbative theoretical predictions PP[. On can hope, however, that the 



five-loop contributions being taken into account will eliminate this inconsistency. 

The authors are grateful to Dr. E. Blagoeva for sending a copy of her important 
work. 
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Table I: Eigenvalue exponents estimates obtained for the Bose (BFP) and the complex 
cubic (CCFP) fixed points at N = 2 and N = 3 within the four-loop approximation in 
e (2e = 1) using Borel transformation with a conformal mapping. 



Type of 
fixed point 


N = 2 


N = 3 


Ai 


A 2 


Ai 


A 2 


BFP 


-0.395(25) 


0.004(5) 


-0.395(25) 


0.004(5) 


CCFP 


-0.392(30) 


-0.030(10) 


-0.400(30) 


-0.015(6) 



Table II: Critical exponents calculated for the Heisenberg (HFP), the Bose (BFP) and 
the complex cubic (CCFP) fixed points at N = 2 and N = 3 within the four-loop 
approximation in e (2e = 1) using Borel transformation with a conformal mapping. 



Type of 
fixed point 


N = 2 


N = 3 


V 


V 


7 


V 


V 


7 


HFP 


0.0343(15) 


0.725(15) 


1.429(20) 


0.0317(10) 


0.775(15) 


1.524(25) 


BFP 


0.0348(10) 


0.664(7) 


1.309(10) 


0.0348(10) 


0.664(7) 


1.309(10) 


CCFP 


0.0343(20) 


0.715(10) 


1.404(25) 


0.0345(15) 


0.702(10) 


1.390(25) 
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